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SUMMARY 


Flat panel flutter at high supersonic Mach number is analyzed using 
NASTRAN Level 16.0 by means of modifications to the code. Two- 
dimensional plate theory and quasi-steady aerodynamic theory are 
employed. The finite element formulation and solution procedure 
are presented. Modifications to the NASTRAN code are discussed. 
Convergence characteristics of the iteration processes are also 
briefly discussed. Effects of aerodynamic damping, boundary support 
condition and applied in-plane loading are included. Comparison of 
nonlinear vibration and linear flutter results with analytical 
solutions demonstrate that excellent accuracy is obtained with 
NASTRAN. 


INTRODUCTION 

Panel flutter is the self -excited oscillation of the external skin of 
a flight vehicle when exposed to an airflow along its surface. The 
classic approach using linear structural theory indicates that there is 
critical (or flutter) dynamic pressure above which the panel motion 
becomes unstable. Since the linear theory does not account for 
structural nonlinearities, it can only determine the flutter boundary 
and can give no information about the flutter oscillation itself. 

A great quantity of literature exists on linear panel flutter (e.g. 
refs. 1 and 2 plus others too numerous to mention). 

For large deflections, the nonlinear effects, mainly due to midplane 
stretching forces, restrain the panel motion to bounded limit cycle 
oscillations with increasing amplitude as dynamic pressure increases. 
Therefore, for realistic assessments and understanding of panel flutter, 
the nonlinear theory should be used. An excellent survey on both 
linear and nonlinear panel flutter through 1970 is given by Dowell 
(ref. 3). 
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To investigate large amplitude panel flutter, a number of approaches 
can be used. A modal approach with direct numerical integration has 
been used by Dowell (refs. 4 and 5). The major disadvantage in using 
this approach is its long computing time. The harmonic balance method 
can be used to determine limit cycles; see for example, Eastep and 
McIntosh (ref. 6) and Kuo et. al. (ref, 7). This approach, however, 
is quite complicated in mathematic manipulations. Morino (refs. 7 and 
8) also used the pertubation method to obtain neighboring solutions to 
the linear problem. 

The finite element method has been used’ successfully in investi- 
gating linear panel flutter (refs. 9 to 15). Because of its versa- 
tile applicability, effects of aerodynamic damping, complex panel 
configuration (e.g. delta planform in ref. 11, and rhombic planform in 
ref. 13), flow angularity, midplane forces, and anisotropic material 
properties can be conveniently included. Recently, the finite 
element method has been applied successfully in large amplitude 
vibrations of beam and plate structures (refs. 16 to 18). Thus, it 
is logical to ^tend the finite element application to study the limit 
cycle oscillations of panels. 

The purpose of this paper is to describe a large deflection super- 
sonic panel flutter capability available for NASTRAN Level 16.0 by 
means of DMAP sequences and modifications of the code. The paper 
Includes a brief discussion of the theoretical formulation and 
solution procedure. Effects of aerodynamic damping, initial in-plane 
loading and boundary support condition are included. DMAP sequences 
required for nonlinear panel flutter analysis and an example of input 
bulk data are given in the Appendices. 
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amplitude of oscillation 
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aerodynamic damping matrix 
modulus of elasticity 
interpolation function 
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SYMBOLS (CONT'D) 


S/ 


1 s -^-1 


[k] 

[k*^] 

[k®l 


M 

[m] 
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N 


xo 



w 


x,y,z 


a 
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SD 


aerodynamic damping parameter, equation (21) 
thickness 

stiffness matrix 
differential stiffness matrix 
geometrical stiffness matrix 
Mach number 
mass matrix 

inplane force due to deflection, tension positive 

applied inplane force, tension positive 

aerodynamic force 
generalized aerodynamic force 
dynamic pressure 
time 

nodal displacements 
flow velocity 
deflection 
coordinates 
damping factor 

norm 

complex eigenvalue, equation (19) 
dynamic pressure parameter 
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SYMBOLS (CONT’D) 


y aerodynamic damping coefficient, equation (8) 


Poisson's ratio 
panel mass density 


Pa 

air mass density 

w 

eigenvector 

= a -f i (1) 

response of system 

U) 

frequency 

0 

reference frequency 

Subscripts: 


aa 

analysis 

ee 

element 


THEORETICAL FORMULATION AND ITS SOLUTION 
Formulation of Matrix Equation of Motion 


The panel is represented by a flat thin plate of unit width in 
bending as shown in figure 1. The transverse dynamic equilibrium 
equation may be written as: 
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is the membrane force Induced by large deflections, and is 
the initial in-plane loading. For sufficiently high supersonic 
speeds (M > 1.6), the aerodynamic pressure can be described by 
the two dimensional aerodynamic theory: 


p(x.y.t) 


2a r 3w ^ 1 M^-2 3w‘ 
3 V M^-1 9t. 


( 3 ) 


In the finite element method, the stiffness equations of motion 
for a plate element under the influence of elastic, initial in— 
plane, large deflection, and inertia forces (ref. 17) with the 
inclusion of aerodynamic forces may be written as: 


([k ] + ] + [k® ]){u } + } = {Q(t)> (4) 

ee ee ee e ee e 

The stiffness [k ], differential stiffness [k*^ ] , and mass 

ee 

fin 1 matrices have been well developed for almost every plate 
ee 

finite element available. The geometrical stiffness matrix 

fk^ 1 has been derived in references 16 and 17 for beam and 
ee 

rectangular plate elements. The development of the aerodynamic 
matrices follows the method proposed by Olson (refs. 9 and 11). 
The virtual work, U, of the aerodynamic force is 


U = Q.u. 

3 3 

= //p(x,y,t) w dx dy 

Assuming the displacements are exponential functions of time 
w(x,t) = w(x)e^^ 


(5) 


( 6 ) 
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where, in general, - a + iw. Substituting expressions for the 
aerodynamic pressure, equation (3), and the displacement func- 
tions, equation (6), the virtual work becomes 



The deflection function for a particular element is usually assumed 
in the form: 


w(x) = Z f, 

j J 


(X) Uj 


{£)’' (u) 


(9) 


where f is the interpolation function corresponding to the 
element nodal degree-of-freedom u.. Introducing the expression 
for w(x) , equation (7) yields ^ 


U =!- 1^ {u^}^ [a^^] {u^} - yft {u^}*^ [d_] {u^}[ 
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where 
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is the non-symmetric aerodynamic matrix and 
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is the aerodynamic damping matrix. The generalized aerodynamic 
forces are 




( 13 ) 
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and their substitution into equation (4) yields the dynamic 
equllibrltim of the panel in the form: 
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The aerodynamic damping matrix, equation (12), can be related to 
the mass matrix by the expression: 


[d 1 * — r [™ 1 
ee ph ee 

and equation (14) takes the final form for a finite element as, 
““ee’ * “‘■'ee' + >V> 

(16) 

+ [m ] + ^ [m ]) (u } - 0 
ee ph ee e 


Solution Procedure 

Assembling the finite elements, applying the kinematic boundary 

conditions, and dividing by.D ^equation (16) leads to a nondlmen- 

^ 3^ 
a 

slonal eigenvalue problem of the form: 
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are the nondlmensional dynamic pressure parameter and eigenvalues, 

respectively. The eigenvalues can be put into more convenient 
form as 
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IS the nondlmensional aerodynamic damping parameter, and 
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is a convenient frequency scale. For typical panels, ranges 
from 0 to 50 approximately, as given in figure 2 of reference 2. 

In determining the eigenvalues K in equation (17) for a given 
dynamic pressure X, the iterative procedure and equivalent 
linearization technique discussed in detail in reference 17 was 
employed. A simple flow diagram of the procedure is shown in 
figure 2. The solution procedure is illustrated briefly as follows. 
For a given X, first the linear flutter problem is solved 

-=I"aal ^ '>‘"aal * ^ '“aa» '*>0 


where represents the linear mode shape normalized by its 
maximum ^ components. The first approximate displacement is 
then expressed in the form 


{u^)^ . c Raal({*}„ e(“ * 


(24) 


where c is a given amplitude of panel oscillations, and a and u) are 
the panel response parameters related to K and g^ by equation (30) . 
An equivalent geometrical stiffness matrix ^^^aa^ 

can be obtained using {u } , and equation (17) is Approxi- 

mated by a linearized ^ ^ eigenvalue equation of the form 


K[m^al ^'*’^1 “ 


] + 

aa 




where K is the eigenvalue associated with amplitude c, and {<()}^ 
is the corresponding mode shape. The iterative process can 
be repeated until a convergence criterion is satisfied as shown in 
figure 2. The maximum displacement norm convergence criterion pro- 
posed in reference 19 was used in the present study and is defined as 
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where Au is the change in displacement component j during iteration 
cycle n, and u is the reference displacement. The reference 

displacement * is the largest displacement component of the 
corresponding "type". For instance in a panel flutter problem in- 
volving deflections and rotations, the reference displacement is the 
largest deflection component and the largest rotation, respectively. 
In addition a frequency norm is also introduced in the present study 
and is defined as 


where Ak^ is the change in eigenvalue during iteration cycle n. 

A typical plot of the maximum and frequency norms versus number 
of iterations for a simply supported panel is shown in figure 3. 

A modified absolute norm and a modified Euclidean norm defined in 
reference 19 were also calculated. They fall in between the maxi- 
mum and frequency norms, and therefore, are not plotted on the 
figure. In the examples presented in the following section, conver- 
gence is considered achieved whenever any one of the norms reaches 
a value of 10“^. 

Equation (17) indicates that when X=0 the problem degenerates into 
large amplitude vibrations of invacuo panels. The matrices [k] , 

tk"*], [k®], and [m] are all symmetric and the eigenvalues are real 
and positive. As X is increased from zero, two of these eigenvalues 
will usually approach each other and coalesce to k at X = X , 
and become complex conjugate pairs 


K 


Kr + i K 


I 


(28) 


for X > X . Here X is considered to be the lowest value of X 
cr cr 

for which coalescence occurs among all limit cycle amplitudes and 
usually corresponds to c = 0. A typical plot of K versus X is shown 
in figure 4. In the absence of aerodynamic damping (g - 0), the 
flutter boundary simply corresponds to X . When X is oelow X , 
any disturbance to the panel decays and (c/h) 0. 

For X > a periodic limit cycle oscillation exists which 

increases in amplitude as X increases. This can be seen more 
clearly by noting that the eigenvalue with a negative Imaginary 
part leads to an instability (see ref. 13) and relating the 
complex eigenvalues to the panel response parameters a and O) 
as follows. Rewrite equations (20) and (28) as 
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(29) 





I + 1 Kj) - 0 


which can be solved for to give 






(30) 



The complete panel behavior is characterized by plotting the varia- 
tion of a + i(jJ with increasing dynamic pressure X. Amplitude 
increases when ct becomes positive. A typical plot is shown in 
figure 5. 


MODIFICATIONS TO THE NASTRAN CODE 

To incorporate this new capability into NASTRAN, four 
existing NASTRAN subroutines must be modified. These subroutines 
are DEAR, REAR, SDRIA (SDRIAZZ on CDC computers because of multiple 
entry points) , and XMPLBD. DEAR was modified in the same way as 
shown in reference 17 . 

Subroutine REAR was modified to calculate the aerodynamic matrix 

[a 1. This matrix is multiplied by the parameter DPMN = 2q/$. 
ee-* 

DPMN is input via a PARAM card in the BULR DATA deck. DPMN is 
passed to REAR through blank common from module EMC. The new EMC 
calling sequence allowing for the DPMN parameter is shown as 
follows: 
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EMG EST,CSTM,MP1;dIT,GE 0M2,/KELM,KDICT,MELM,MDICT, ,/V,N,NOKGGX/ V, 
N,NOMGG/C,N,/C,N,/C,N,/C,Y,COUPMASS/C,Y,CPBAR/C,Y,CPROD/C,Y, 
CPQUADl/ C , Y , CPQUAD2/C , Y , CPTRIAl/C , Y , CPTRIA2/ C , Y , CPTUBE/C ,Y , 
CPQDPLT/C , Y , CPTRPLT/C , Y , CPTRBSC/V, Y , DPMN $ 


The default value for DPMN, which is set in XMPLBD, is zero (0) , 
This means if the PARAM card for DPMN is omitted, [a ] will 
make no contribution in equation (14). 


Subroutine SDRIA was modified to calculate the real part of 

{^}n complex eigenvector generated b> 

the module CEAD. To avoid entering the modified section of code 
each time SDRIA is called, a new parameter, IFLUT, was added to 
the DMAP calling sequence for module SDRl. The contents of IFLUT 
are passed through blank common from SDRl to SDRIA. The default 
value for IFLUT, which is set in XMPLBD, is zero (0). When IFLUT = 0, 
the new code in SDRIA will not be executed. To set IFLUT = 1 and 
execute the new code in SDRIA, the following calling sequence for 
the SDRl module is used: 


SDRl USET,PHIA, ,,G0,GM,KFS/PHIG,, BQG/1/*REIGV1 $ 

The underlined parameter sets IFLUT to 1. 

Once the changes were made to DBAR, KBAR, SDRIA, and XMPLBD, they 
were compiled and replaced the old DBAR, KBAR, SDRIA, and XMPLBD 
in the NAS TRAN object library. Link 1, Link 3, Link 12, and Link 13 
were relinked, creating a new executable NASTRAN. Although this 
procedure was done on a GDC computer, similar procedures will produce 
similar results on the IBM and UNIVAC computers. 

To use this capability in NASTRAN, the DMAP sequence shown 
in Appendix A must be used. This sequence uses many of the new DMAP 
convenience features in Level 16 of NASTRAN. One of the features 
allows the KEPT module to have a variable parameter. The variable 
parameter NL00P is used for KEPT in this DMAP sequence. NL00P is 
input on a PARAM card in the BULK DATA deck. It sets the maximum number of 
iterations of the inner loop shown in figure 2. The only other input 
required to use this capability is the addition of another PARAM card 
in the BULK DATA deck. The parameter AMP corresponds to c and is 
used to specify the amplitude of vibration of this structure. This 
capability was added to an in-house version only and is not 
available in any standard NASTRAN level. 
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RESULTS AND DISCUSSION 


The large deflection panel flutter analysis developed for use with 
NASTRAN has been applied to various panels. A typical BULK DATA 
deck for a simply supported panel at — = 0.6 and X = 600.0 is 
given in Appendix B. 


Convergence Study 

Numerical results for the first two eigenvalues at X = 0 and 
for the coalescence for a simply supported panel and a clamped 
panel are shown in Table 1. The exact results for eigenvalue 
coalescence are from reference 20. It is seen that an excellent 
approximation to the exact results is obtained with only eight 
elements . 


The influence of large deflections on in-vacuo frequencies for a 
simply supported panel is given in Table 2. Analytical solutions 
using three different approaches from reference 21 are also given. 
Comparison of the NASTRAN results with the reference 21 methods show 
that the eight-element approximation gives very good results. There- 
fore, eight elements were used in modeling the panels in all the 
flutter results presented. 

Simply Supported Panel and Effect of Aerodynamic Damping 

Plots of the eigenvalues verses dynamic pressure for a simply supported 
panel at two different panel amplitudes ^ “ ^.0 (linear theory) and 
0.6, are shown in figure 4. The complete panel behavior is charac- 
terized by plotting the (ot + ico) variation with increasing dynamic 
pressure X, using equation (30) and figure 4, as shown in figure 5. 

For the case of negligible aerodynamic damping, instability 

does not set in until after the two undamped natural frequencies have 
merged. If some damping is present, the instability sets in at a 
somewhat higher value as indicated in figure 5. This occurs when 
a = 0 in equation (30). By routine algebraic manipulation, this 
instability occurs at the value of X when 

(32) 



and the corresponding limit cycle frequency is 
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However, as discussed earlier, this instability is not catastrophic* 
The panel response does not grow indefinitely, but rather a limit 
cycle oscillation is developed with increasing amplitude as X 
increases . 


Boundary Support Effect 

In figure 6, the panel amplitude of the limit cycle oscillation is given 
as a function of X for various panel edge restraints. The most interest- 
ing result is that the limit cycle motions are different for hinged- 
clamped and clamped-hinged panels. This occurs because the aero- 
dynamic matrices are different for the two support conditions, 
which leads to different deflection shapes for the panels as well 
as different geometrical stiffness matrices. 


Effect of In-Plane Loading 

Panel amplitude versus X for severa] applied in-plane forces acting on a 
simply supported panel is shown in figure 7. The classical Euler 
buckling load for simply supported panels is N^r “ -TT^D/a^. The 
total membrane force is composed of the applied in-plane load 
and the membrane force N induced by large deflections of the 
panel.- Figure 7 shows that the applied compressive in-plane force 
reduces the critical dynamic pressure. However, as the dynamic pressure 
is increased the panel amplitude increases, which induces tensile in- 
plane forces that counteract the applied compressive forces. This 
process continues until a flutter dynamic pressure is reached which 
corresponds to a given limit cycle amplitude. 


CONCLUDING REMARKS 

A large amplitude supersonic panel flutter capability has been 
developed for use with NASTRAN Level 16.0 by means of DMAP 
sequences and modifications to the code. An aerodynamic matrix 
for a two-dimensional plate element has been developed for NASTRAN 
by modifying subroutine KBAR. The iteration process has been 
implemented in NASTRAN through PARAM NL00P in bulk data deck, 
modifications in subroutines DBAR, KBAR, and SDRIAZZ, and the DMAP 
sequences. Examples which include effects of aerodynamic damping, 
applied inplane forces and various support conditions have demon- 
strated the versatility of the method. 
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APPENDIX A 


DMAP SEQUENCES 


ID 

APP 

BEGIN 

XDMAP 

FILE 

GPl 

SAVE 

CHKPNT 

GP2 

CHKPNT 

PARAMU 

PURGE 

COND 

PLTSET 

SAVE 

PRTMSG 

PARAW 

PARAM 

COND 

PLOT 

SAVE 

PRTMSG 

label 

CHKPNT 

GP3 

CHKPNT 
TA 1 

SAVE 

COND 

PURGE 

CHKPNT 

PARAM 

PARAM 

EMG 


SAVE 

CHKPNT 

COND 


nlpf^twod 

dmap 

$ 

G0«ERR = 2 tL I ST • 

LAMA=APPEND/PH I A=APPEND $ 

GEOMl « GEOM2 « /GPL » EQEX IN » GPDT « CSTM • bGPDT ♦ S 1 L/V tN t LUSET/ V *N» 
NOGPDT $ 

LUSET S 

GPL f EQEX I N ♦ GPDT ♦ CSTM , BGPDT f b I L S 
GE0M2%EQEXIN/ECT $ 

ECT $ 

PCDB//C »N ♦PRES/C #N« /C *N » /C* N* /V •N^NOPCDB S 
PLTSETXf PLTPAR,GPSETS«ELSETS/N0PCDB S 
Pl*NOPCDB S 

PCDB tEOEX IN *£01 /PL I ^dL I X » PL i PAR ,GP:r3L \ t EL::>b i^/v«N^N:3lL/ VtN* 
JUMPPL0T = “1 5> 

NSiLf JUMPPLOT S 

pltsetx// S 

//C*N*MPY/V»N,PLTFLG/C#N« 1/C *N# 1 S 
//C ♦ N f MPY/V ♦ N • PF I LE/C f N 1 0/C ♦ N ♦ 0 S 
PlfJUMPPLOT 

PLTPAR fOPSETS^ELSETS^CASECC^ BGPDT •EQEX IN fSiL# « • •/PlOTxI/ V«N* 
NSlL/VfNfLUSET/VWJ# JUMPPLOT/V •Nf PLTFLG/V*N«PF ILE S 
JUMPPLOT • PLTFLG • PF I LE $ 

PLOTXI// S 
PI s 

PLTPARf GPSETS^ELSETS S 

GE0M3«EQEX IN • GE0M2/SLT • GPTT/V • N • NOGRA V S 
SLTtGPTT S 

ECT • EPT ♦ BGPDT IL fGPTT • CS I M/ tib 1 • Gt I • GPEC 1 •/V • Nf LUbE I / V»N» 
NOSIMP/C»N» 1/V*N*N0GENL/V#N»GENEL 

nosimp*nogenl*genel S 
error I •NOS I MP $ 

OGPST/GENEL 4 

ESTfGPECTfGEI .OGPST S 

//C , N ♦ ADD/V • N ♦ NOKGGX/C , N , 1 /C • N , 0 4 

//C • N • ADD/V • N , NOMGG/C • N • 1 /C • N • 0 4 

EST»CSTM,MPTfDlTf GE0M2 « /KELM »KDICT •MELM^MDICT* • /V »N fNOKGGX/ V « 

N* NOMGG/C •Nt /C »N • /C • N, /C • Y • CO JPMASS/C • Y , CPBAR/C « Y •CPROO/C • Yf 

CPQuADl/CfYf CPOUAD2/C, YfCPTR I Al/C*YtCPTRl A2/ C , Y • CP I uBE/C , Y , 

CPQDPLT/C • Y • CPTRPLT/C • Y • CP 1 RBbC/V • Y • DPMN S 

NOKGGXfNOMGG 4 

KELMfKDICT^MELMfMDl CT S 

JMPKGGfNOKGGX 4 
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EMA 

CHKPNT 

LABEL 

COND 

EMA 

CHKPNT 

COND 

GPWG 

OFP 

LABEL 

EQUI V 

CHKPNT 

COND 

SMA3 

CHKPNT 

LABEL 

PARAM 

GP4 


SAVE 

COND 

PARAM 

PURGE 

CHKPNT 

COND 

JUMP 

LABEL 

COND 

GPSP 

SAVE 

COND 

OFP 

LABEL 

EQUIV 

CHKPNT 

COND 

MCEl 

CHKPNT 

MCEE 

CHKPNT 

LABEL 

EQUIV 


GPECT «KD I CT«KELM/kGGX * GPST S 
KGGX.GPST S 
JMPKGG S 
ERRORS. NOMGG $ 

I CT . MELM/mGG ./C.N. — 1/C«Y.WTMASS=1 «0 $ 

MGG S 

LBLl .GRDPNT S 

BGPDT • CSTM ♦ t-QhX i N • Mbb/OGPW^ / V ♦ Y ♦ bKUHN 1 / C * Y ♦ W I 
OGPWG « • « « f // s 
LBLl S 

KGGX •KGG/NOGENL $ 

KGG S 

LBLllfNOGENL $ 

GEI ♦ KGGX/KGG/V f N «LUbE f / V ♦ N t NOGENL / V • N « NOb I MP S 
KGG $ 

LBLl 1 $ 


//C.N.MPY/V.N.NSK IP/C.N.O/C.N.O $ 

CASECC,GE0M4,EQExIN,SIL.GPD 1 .BGPDI .C-SIM/RG. Y:.,UiEl .Ai>E l/V.N. 
KcJtl l/V.N.REACT/V.N. 

NSKIP/v.N,REPEAT/V.N,NOSET/V,N,NOL/V.N,NOA/C,y.SUB1D 1 

^ ' NGLE , OM I T , RE AC T , N£K I P , REPE A T , NOSE T , NOL . NOA 4 

cWkOko « NOL S 

//C ♦ N • AND/V/ * N* NOSR/V « N ♦ S I NGLE/V • N ♦ RE AC I S 
GM.RG.GO.KOO.LOO.PO.UOOV.RUOV, YS.PS.KFS.KSS.USET.ASET.QG S 

LBL4D. REACT 4 
ERRORS S 
LBL4D S 
LBL4.GENEL t 

GPL.GPST.USET . S 1 L/OGPST/V «N « NOGPST S 

NOGPST S 

LBL4. NOGPST Sb 

OGPST,,,,,// 4 

LBL4 4 

KGG«KNN/MPCF1 /MGG.MNN/MPCFI 

KNN.MNN S 
LBL2.MPCF2 S 
USET.RG/GM S 
GM 4 


USET.GM.KGG.MGG. . /KNN.MNN.. S 
KNN.MNN S 
LBL2 4 

KNN.KFF/SINGLE/MNn.MFF/SINGLE 4 
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chkpnt 

COND 

SCEl 

CHKPNT 

label 

EQUIV 

CHKPNT 

COND 

SMPl 

CHKPNT 

SMP2 

CHKPNT 

LABEL 

RBMG2 

CHKPNT 

SSG 1 

CHKPNT 

EQUIV 

CHKPNT 

COND 

SSG2 

CHKPNT 

LABEL 

SSG3 

SAVE 

CHKPNT 

COND 

MATGPR 

MATGPR 

LABEL 

SDRl 


CHKPNT 

SDR2 

PARAM 

OFP 

SAVE 

COND 

PLOT 

SAVE 


KFFfMFP S 
LBL3tS INGLE ^ 

USET •KNNfMNNt t/KFE ♦KFS«KSS«MFF*f % 

KFS«KSS*KFF f MFF S 
LBL3 S 

KFFfKAA/OMIT/MFFfMAA/OMIT ^ 

KAA«MAA S 
LBLBfOMIT S 

USET «KFF* • ♦/GOf KAA«KOO»LOO«UOO* t • ♦ ^ 

GOf KAAfKOOfLOOtUOO 5» 

uset#go#mff/maa S 

MAA S 

LBL5 S 

KAA/LLL S 

LLL S 

SLT • BGPOT *CSTM«SIL* tST ♦ MPT f GPTT ♦ EDT < MGG • CAStCC « D I T /PG/ V «N ♦ 
LUSET/C^Nfl h 
PG S 

PG# PL/NOSET % 

PL S 

LBL10#NOSET S 

USE T#GM#YS#KFS#G0## PG/ f PO # PS • PL S 
PO#PS#PL S 
LBL 1 0 S 

LLL# KAA^PL #LOO#KOO *PO/ULV #UOOV *RULV ♦ RUOV/V «N»OMlT/V#Y# 1RES = -1/ 

C • N • I /V # N # EPS I S 
EPS I S 

UL V • UOOV f RULV « RUOV S 
LBL9#IRES S 

GPL#USETfSIL#RULV//C#N#L S 
GPL#USETfSIL#RUOV//CfN#0 S 
LBL9 S 

USET f PG# ULV # UOOV • YS # GO # GM • PS • KFS #KSS # /UGV ♦ PGG % QG/C «N* 1/C«N, 

BKLO $ 

UGV#QG#PGG S 

CASECC fCSTM#MPT#DlT#EQEXlN#i5l L #GPTT #EDTf6GPDT## QG#UGV # EST ♦ »PGG/ 
OPGl #OQGl #OUGVl #0ES1 #OEFl fPUGV I /C #N#BKLO S 
//C # N # MP Y/V # N ♦ CARDNO/C # N « 0/C # N • 0 $ 

OUGV 1 # OPG I # OQG I # OEF 1 • OE S 1 # //V # N # C ARDNO S 
CARDNO S 
P2#JUMPPL0T 3, 

PL T PAR « GPSETS # ELSET S • CASECC # BGPDT #t QEX I N • S I L # PUGV 1 # • GPECT # OE SI/ 
PLOTX2/V # N # NS I L/V #N tLUSE T/V t N « JUMPPLOT / V # N ♦ PLTFLG/V # N # PF 1 LE S 

PFILE S 
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PRTMSG 

UABEL 

TAl 

OSMGl 

SAVE 

CHKPNT 

DPO 


SAVE 

COND 

CHKPNT 

PARAM 

LABEL 

ADD 

CHKPNT 

CEAD 

OFP 

COND 

SDRl 

ADD 

DSMGl 

CHKPNT 

ADDS 

CHKPNT 

EOUI V 

CHKPNT 

COND 

MCE2 

CHKPNT 

LABEL 

FQUIV 

CHKPNT 

COND 

SCEl 

CHKPNT 

label 

EQU 1 V 

CHKPNT 

COND 

SMPl 

SMP2 


PLOTX2// S 
P2 t 

ECTfEPT *BGPDT ♦ S I L » GPDT ♦ CSTM/X 1 ,X 2 
NOS I MP/C « N« 0 /V *N»N 0 GENL/V »N« GENEL 
CASECC * GPTT ♦ S 1 L ♦ EDT « UGV « CSTM « MPT « 
DSCOSET S 
DSCOSET S 
KDGG $ 


♦ECPT«GPCT/V«N*LUSET/ VtN« 

S 

ECPT ♦ GPCT * D 1 T/KDGG/ V«N» 


DYNAMICS, GPL, EEO , EQDYN/ V , N , 

LUSET/V,N,LUSETD/V,N,N0TFL/V,N,N0DLT/V,N, NOPSDL/V,N,NOFRl/ 
n»nonlft/v,n,notrl/v,n,noeed/c ,n,/v,n,noue $ 

NOEED S 

errors, NOEED £ 

EEC $ 


//C ,N , MPY/V,N, NE I GV/C , N , 1 XC , N , - 1 $ 

NLVIB $ 

KAA ,/KTT/C ,N , (-1,0,0«0)/C,N,<0.0,0»0) $ 
KTT $ 


V, 


KTT, *maa,eed,casecc/phi A,LAMA,0EIGS/S,N,NEIGV $ 

OEl GS ,LAMA//S , N, CARDNO £ 

ERR0R4,NEIGV £ 

USET, ,PHIA, , ,GO,GM, ,KFS/PH1G, , BOG/ 1 /*RE T G*/ 1 £ 

PH I G/PHl AMP/V f Y ♦ AMP S 

CASECC, ,S1L, *PH1 AMP,CSTM,MPT,ECPT,GPCT,0I T/KNGG/DSCOSET/1 £ 

KNGG S 

KGG , KDGG , KNGG , , /KSGG £ 

KSGG £ 


KSGG , KSNN/MPCF2/MGG , MSNN/MPCF 2 £ 
KSNN,MSNN £ 

LBL23,MPCF2 £ 

USET, GM, KSGG, MGG/KSNN,MSNN £ 
KSNN,MSNN £ 

LBL2S £ 


KSNN ,KSFF /S I NGLE/MSNN , MSFF/S I NGLE £ 
KSFF,MSFF £ 

LBL3S, SINGLE £ 

USET ,KSNN,MSNN/KSFF,KSFS, ,MSFF £ 
KSFF,KSFS,MSFF £ 

LBL3S £ 

KSFF,KSAA/OMlT / MSFF , MSA A/OM I T £ 
KSAA,MSAA £ 

LBL5S,0M1T $ 

USET ,KSFF/GS0,KSAA,KS00,LS00,US00 £ 
USET ,GSO, MSFF/MSAA £ 
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chkpnt 

LAL3EL 

COPY 

COPY 

REPT 

SDR2 

OFP 

SAVE 

COND 

PLOT 


SAVE 

PRTMSG 

LABEL 

JUMP 

LABEL 

PRTPARM 

LABEL 

PRTPARM 

LABEL 

PRTPARM 

LABEL 

PRTPARM 

LABEL 

PRTPARM 

LABEL 

PRTPARM 

LABEL 

END 


k*saa »msaa s 

LBL5S S 

KSAA/KAA/ I PARM=- 1 s 
MSAA/MAA/ JPARM=- 1 * 

NLV IB • NLOOP S 

CASECC ♦ CSTMfMPT*D I T ^ EQEX IN*S IL f « *3GPDT , L AM A ♦ BOG ♦ PH I G t E ST t #/• 

OBQGl fOPHlGiOBESl ,OBEFl , PPH I G/C f N ♦ RE I G S 

CPHIGfCBQGl ,CBEF1 *CBES1 ^ 1//V.NSCARDNO S 

CARDNO S 

P3 • JUMPPLOT S 

PLTPAR»GPSETSfELSETS«CASECC*BGPDT «EQEXlN*SIL« fPPHlG«GPECT « 

ORESl/PLCTX3/VfNfNS IL/V*NfLUSET/V*N, JUMPPLOT/V«NfPLTFLG/V*N, 

PFILE S 
PFILE S 
PL0TX3// 3> 

P3 £ 

FINIS S 

error 1 S 

//C«N»-1/C#N»NMDS % 

ERROR2 $ 

//C •Nf~2/C *N»NMDS $ 

ERR0R3 S 

//C«Nt-3/CfNfNMDS S 
ERR0W4 $ 

//CfN«-4/C*NfNMDS S 

errors $ 

//C»Nf-~5/C%NfNMDS $ 

ERRORS $ 

//C»Nf-6/C*NfNMDS S 
FINIS S 
S 
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APPENDIX B 

INPUT BULK DATA CARDS 


S 

s geometry 

GROSET 
GRID 1 

GRID 2 

GRID 3 

GRID 4 

GRID 5 

GRID 6 

GRID 7 

GRID 8 

GRID 9 

GRID 20 

SPC 1 

S 

s STRUCTURAL 
BAROR 


AND CONSTRAINTS 


0*0 
0*125 
0*25 
0.375 
0.50 
0.625 
0.75 
0.875 
1 .0 
0.0 
13 


246 


10. O 

0.0 


123456 

0#0 


AND AERODYNAMIC 
15 


ELEMENTS 


CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
MATl* 

*MT1 
PARAM 
S DPMN 
S 
S 

PARAM 
PBAR* 

*PB1 
S 

s CONTROL 


1 I 

2 2 

3 3 

4 4 

5 5 

6 6 

7 7 

e 8 

25 1 .0 

0.4367901341 
COUPMASSl 
= 2.0J^Q/ (BETA ) 

WHERE Q * RH0*V**2/2.0. 
BETA = SORT (MACH 


20 


2 

3 

4 

5 

6 

7 

8 
9 


» dynamic 

N0.**2 - 


DPMN 
15 
1 .0 


DATA 


EIGC 
+ INVl 
S AMP 
S 
S 

PARAM 

PARAM 

S 


1 

41, 


600*0 


INV 

-12* 


25 


MAX 
41 .0 


PRESSURE 

1 .0 ) 

2.289429 


1 #0 


-14.0 1.00 

AMPLITUDE/SQRT( I/A) = SORT ( 1 2 . 0 ) »C/H 
WHERE I = AREA MOMENT OF INERTIA 
A = AREA 

AMP 2.0784610.0 

NLOOP 3 


1 1 
SQRT( 12.0 )*0»6 


2.07846 


* applied inplane 

FORCE* 1 

+FCE -1.0 0, 

ENDDATA 


loading 

9 

0 0.0 


9.86960440109 


♦ MTl 


♦PBl 


+ INV1 


+FCE 
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TABLE 1. IN VACUO EIGENVALUES AND COALESCENCE RESULTS 
FOR SIMPLY SUPPORTED AND CLAMPED PANELS 


S^mplv Supported Panel 
In Vacuo Coalescence 


Number of 
Elements 

^1 

^2 

X 

cr 

K 

cr 

2 

98.1795 

1920.00 

398.536 

1206.32 

4 

97.4597 

1570.87 

342.347 

1043.47 

8 

97.4123 

1559.35 

343.280 

1051.22 

Exact 

97.4091 

1558.55 

343.3564 

1051.797 

(ref. 20) 







Clamped Panel 



In 

Vacuo 

Coalescence 

Number of ^ 

Elements 1 

•^2 

X 

cr 

K 

cr 

2 

516.923 

6720.00 

922.388 

3618.46 

4 

501.894 

3874.23 

636.437 

2721.38 

8 

500.648 

3808.34 

636.586 

2740.16 

Exact 

500.564 

3803.54 

636.5691 

2741.360 


(ref. 20) 
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for each element 


assemble 


SET j - 1 


LOOP ON A. 
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SET ITERATION COUNT n = 1 
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COMPUTE [k® 1 

aa eg 
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CONVERGENCE 

TEST 
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DAMPING FACTOR 








97 


Figure 7. Limit cycle amplitude versus dynamic pressure for simply supper 
different in-plane forces (g = 0). 



